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Abstract 

In this paper we introduce a model describing diffusion of species by a suitable 
regularization of a “forward-backward” parabolic equation. In particular, we prove 
existence and uniqueness of solutions, as well as continuous dependence on data, for 
a system of partial differential equations and inclusion, which may be interpreted, 
e.g., as evolving equation for physical quantities such as concentration and chemical 
potential. The model deals with a constant mobility and it is recovered from a 
possibly non-convex free-energy density. In particular, we render a general viscous 
regularization via a maximal monotone graph acting on the time derivative of the 
concentration and presenting a strong coerciveness property. 
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1 Introduction 


The model we are introducing may be applied to different situations, dealing with dif¬ 
fusion of different species (located in some domain C R"' and) described in terms of 
concentration. In the following, moving from the classical approach in thermodynamics 
which leads to the well-known Cahn-Hilliard equation, we introduce our point of view 
and make some comments on its thermodynamical consistency. In particular, we are fo¬ 
cusing on diffusion in solids and we have in mind, as a possible hnal application, hydrogen 
storage in metals. 
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C.N.R. Pavia. GT acknowledges partial financial support from the Italian GNFN (Gruppo Nazionale per 
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Regularization of a forward-backward parabolic equation 


Models for species diffusion. In the simplest setting, a mathematical model describing 
species diffusion in a solid can be obtained by combining the following three ingredients: 

1) the mass-balance law for the concentration u 

h + divh = 0; (1.1) 

2) a linear constitutive prescription relating the flux of diffusant h to the gradient of 
chemical-potential /i through a mobility constant m: 

h=— mV/i, m > 0; (1.2) 

3) a possibly nonlinear relation between pi and u\ 

fi = 'if\u), (1.3) 

dictated by the derivative fj' of a coarse-grain free energy fj. 

When diffusion is accompanied by phase separation, in general if is assumed to be a 
non-convex function, a typical choice being the double-well polynomial potential: 

if{u) = k{u — /c > 0. (1.4) 

This kind of specihcations renders the system fll.ll) - fll.3p backward parabolic (e.g., for 
01.41) . in the region where 0 < m < 1), an undesirable feature which calls for a suitable 
regularization. 

In order to make the system well-posed, the most popular approach is the so-called 
elliptic regularization of 01.31) : 


p = fjfu) — aAu, a > 0, (1.5) 

which results into the celebrated Cahn-Hilliard system [S]. Yet, other choices have been 
considered: Novick-Cohen and Pego [IB] and Plotnikov [121 have investigated the viscous 
regularization: 

/i = fifu) -|- ail, a > 0; (1.6) 

furthermore, a combination of energetic and viscous regularization, namely, 

/i = fifu) — aAu -\- ail, 

which leads to the so-called viscous Cahn-Hilliard eguation, was derived by Novick-Cohen 
in im and analytically investigated by Elliott and Garcke in na and by Elliott and 
Stuart in [T3|, with computations being carried out in [2]. Its vanishing-viscosity limit 
was studied in [23] • More sophisticated generalizations of the Cahn-Hilliard system that 
still incorporate a viscous contribution have been proposed and investigated in [T61I22] . 
In the same spirit, a similar viscous regularization has been introduced in the paper [1], 
which deals with phase separation in binary alloys driven by mechanical effects. 

The elliptic regularization as a microforce balance. Derivations of diffusion models 
a la Cahn-Hilliard have been proposed by Curtin in [15] and by Podio-Guidugli in |2nj . 
Central to these derivations are the following ingredients: 
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• a system of microforces, distinct from Newtonian forces, which obey their own bal¬ 
ance laws and whose power expenditure is associated to the evolution of the observ¬ 
able fields of interest — in the present case, the concentration m; 

• a collection of constitutive prescriptions, which relate microforces to the actual 
evolution of the system; 

• a dissipation principle, which sifts away thermodynamically inconsistent constitutive 
prescriptions. 

Within this common framework, the elliptic regularization fll.Sp is a particular instance 
of the balance statement: 

div^-|-TT -|-7 = 0, (1.7) 

an instance that arises on adopting the constitutive prescriptions: 

^ = aVu, = ^ — 7 = 0 , ( 1 . 8 ) 

for the vectorial microstress $, and the two scalar-valued fields, tt and 7 , respectively, the 
internal and the external microforce. 

Other constitituve prescriptions may be taken into consideration of course, provided 
that they are thermodynamically consistent. In this respect, two distinct options are 
offered in na and m- In this paper we opt for the former, where thermodynamical 
consistency is embodied by the dissipation inequality: 

0 < (p — 7r)u + ^ ■ Vu — h ■ Vp, (1.9) 

with 0 the free-energy density. 

A standard argument [lO] shows that consistency with fll.9p rules out any constitutive 
dependence of free energy on the time derivative of m; accordingly, one assumes that the 
free energy and the concentration fields are related through a prescription of the form: 

0 = 0(m,Vm). (1-10) 

A consequence of fll.9p and fll.lUl) is that, on introducing the equilibrium parts 

r'* := ^^(w, Vm) and := p - ^{u,Vu), ( 1 . 11 ) 

of, respectively, microstress and internal microforce, the test for consistency of a certain 
constitutive choice with the dissipation inequality fll.9p boils down to verifying that, for 
whatever process, the non-equilibrium parts 

^ ;= tt - 71"*^ (1.12) 

of microstress and internal microforce satisfy, together with the flux of diffusant h, the 
redueed dissipation inequality. 


0 < • Vh - h • V/i. 


(1.13) 
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Regularization of a forward-backward parabolic equation 


In particular, the constitutive prescriptions fll.Sp follow from fll.lip - fll.12p on taking 


Vu) = -a| Vup + ^(u), 


and on choosing 
a choice consistent with fll.131) . 


0 , 


TT 


0 , 


(1.14) 


The viscous and the “non smooth” regularizations. The viscous variant fll.5l) is 

arrived at in a similar fashion: first, we exclude microscopic contact interactions from 
the picture by letting ^ = 0 , and we set to null the external microforce 7 , so that the 
microforce balance fll.7l) reduces to 

TT = 0. (1.15) 

Then, consistent with this choice, we rule out the dependence of free energy on concen¬ 
tration gradient by letting 

0(m, Vu) =(1-16) 
so that the second of fll.lip specializes to 

TT®'^ =/i —(1-17) 


As to the “non-equilibrium part”, the simplest constitutive choice is 

7 r°® = —ail, a > 0, (1-18) 

so that, bearing in mind the second of fll.121) . we recover fll.5p from fll.151) and fll.lTp . 

Note that, actually, the above relation could be introduced in terms of a “dissipation 
functional” <h: the so-called pseudo-potential of dissipation by Moreau, which is a non¬ 
negative and convex functional, equal to zero for null dissipation. In particular, letting 
<h(-u) = f Ihp, the non-equilibrium part could be introduced as 


TT = — 


du 


Note in particular that the assumption on $ and fll.lbp (see fll.lSp l lead to 

7r^% < 0 , 


(1.19) 


( 1 . 20 ) 


which is important to ensure thermodynamical consistency. Indeed, if fll.20p and fll.2p 
hold, then the reduced dissipation inequality fll.ldp is satisfied (recall that = 0 , which 
follows from fll.lip i . fll.l2P i . and fll.lbp ). 

Moreover, we could generalize fll.lSp by a suitable choice of $. Following a suggestion 
in [ 21 ], in this paper we replace fll.lSp with the following prescription; 


(tt”® -|- ail) e /3(ii), 


( 1 . 21 ) 


where, as before, a > 0 is a constant and /3 : M M is a set-valued mapping whose graph 
is (maximal) monotone and containins the origin: 


0 6 /3(0), 


(1.22) 
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as in the case (3 = d( is the sub differential of a non-negative convex function (, with 
C(0) = 0, 

As to the assumption on /3, it is easily seen that the constitutive prescription fll.211) is 
consistent with the dissipation inequality fll.ldp : indeed, owing to the monotonicity of /9, 
we have 

vi e /3{wi) &: V 2 e (3 {w2) ^ (fi - V 2 ){wi - W 2 ) > 0; (1-23) 

thus, if the pair (71°®, it) is compliant with then we can take Vi = —(71°® -|- ait) and 

tci = ft as tests in fll.23p : meanwhile, fll.22p entitles us to choose W2 = V2 = 0; with these 
choices, fll.2311 yields — (tt”® -|- au)u > 0, which entails fll.201) . 

The system we investigate. In order to assemble the system we study, it remains for 
us to: 

1 ) rewrite the balance equation fll.ll) in the light of the constitutive prescription fll.21) 
for the flux of diffusant, whence the partial differential equation: 

u = mA/i. (1-24) 

2) combine the microforce balance, in its form 01.151) . with the second of 01.12p and 
with the prescriptions 01.17P and 01.211) . whence the pointwise inclusion: 

jj, G+ ail +/3{u), (1.25) 

where, with slight abuse of notation, we write p + (3{q) to denote the set {r G M : r — p G 

m}- 

We finally complete the system 01.24l) - 01.25p with a prescription of the initial con¬ 
centration m( 0) (see 01.37cl) below) and a time-dependent prescription of the chemical 
potential on the boundary: 


/i(-,f) = p\){-,t) on r := dfl. (1.26) 

Dirichlet-type boundary data involving chemical potential, rather than the flux of diffu¬ 
sant, are typical in applications. An example in the context of the mathematical modeling 
of hydrogen storage [SHZIEI] is the following: consider that the region where diffusion 
takes place represents a body immersed in a gaseous reservoir at uniform pressure Pg and 
temperature Tg, both of which may be possibly time dependent. The chemical potential 
of the diffusant in the reservoir is related to its pressure and temperature through the 
formula Pg = po-\-RTglog{pg/po) (see for instance [H Eq. 5.16]), where po is the standard 
chemical potential, R is the Avogadro constant, and po is the standard pressure (typically, 
Po = Ibar). If local equilibrium prevails, p is continuous across the boundary E, and hence 
fll.26p holds with p\,{x,t) = Pg{t). Another example is provided by mechanical theories 
that describe the behaviour of a permeable elastic solid immersed in an incompressible 
fluid (see for instance m) ; in these theories, the chemical potential of the fluid is given 
hy pf = po-\- ^{pf — Po), where n is the molar volume of the solvent and p/ is its pressure. 


Energy and dissipation. The precise assumptions we make on the coarse-grain free 
energy -0 are stated in fl2.ll) below. Note in particular that we can allow ijj to be nonconvex 
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Regularization of a forward-backward parabolic equation 


(a feature that, as already pointed out, allows for phase separation) and that we can 
include logarithmic type potentials. On the other hand, we are not able to deal, e.g., with 
sub differential of indicator functions of closed intervals. 

As to the mapping /3, whose choice together with that of the constant a affects the 
dissipative structure of the system, we assume in fll.7p below that it is the subdifferential 
of a non-negative, convex, lower semicontinous potential C (without setting any restriction 
on the growth of Q with (^(0) = 0 . Besides the trivial case /3 = 0, which leads to the 
PDE considered in [18], other possible choices are: 


(i) /3(r) = /So sign(r), where /So > 0 and the sign graph is dehned by 


r {-hi} if r > 0, 

sign(r) = < [-1, -hi] if r = 0, 
[{—1} if r < 0, 


(1.27) 


which, as we shall discuss below, may induce hysteresis. 

(ii) the sub differential of the indicator function I[a,b] of a closed bounded interval [a, b] 
(in our computations we require 0 G [a, b]) 


fjO} if re {a,b), 

/3(r) = dl[a,b]ir) = I [0, -hoo) if r = b, 

[ (—oo, 0] if r = a, 


(1.28) 


forcing the rate of change of concentration to be bounded in the interval [a, b], 
(hi) the subdifferential of the indicator function of [0, +cxo), namely. 


/3(r) = a/[o,+oo)(L) 


{0} if r G (0, -hoo), 

(—cx), 0] if r = 0, 


(1.29) 


which is a choice particularly interesting, for it entails that the concentration at a 
given point cannot decrease, that is, ii in fll.25p has to remain non-negative. 


Remarks on hysteresis. For the viscous variant of the Cahn-Hilliard system — namely, 
the system that arises from fll.61) — it is known that hysteresis (in the sense of irreversibil¬ 
ity emerges in the vanishing-viscosity limit: 

cn —^ 0 , 

provided that is a. non-convex function [HISS]. It is not hard to construct examples 
showing that our constitutive assumptions lead to hysteresis as well, even if ip is convex. 
In order to provide an illustration, we choose fd as in (i) above and we take ip{u) = \ku^, 
with /c > 0. As a result, the system fll.24p - fll.25p becomes: 


dtu = mAfi, 

p G a dtU -h (do sign {dtu) -h ku, 


(1.30a) 

(1.30b) 
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where dt denotes the (partial) derivative with respect to t. We supplement fll.SOap with 
a boundary condition of the form 

= onh, (1.31) 

and we investigate the formal limit when the characteristic time r tends to inhnity, which 
corresponds to the regime of a slowly-varying chemical potential imposed at the boundary. 

On replacing t with the dimensionless variable 


s : = 


t 

r 


and on considering that the (3 in (i) is invariant under time reparametrization, we can 
rewrite fll.dUp as 


r ^dgU = mA/i, (1.32a) 

/i G T~^a dsU + /So sign {dgu) + ku. (1.32b) 

Formally, as r ^ +cxo, the parabolic system fll.321) degenerates into the elliptic system: 

0 = mA/i, (1.33a) 

/i G /So sign ((9sii) + ku. (1.33b) 

Having stipulated with fll.311) that p is spatially constant on the boundary, the homoge¬ 
neous elliptic equation fll.33ap entails that p is spatially uniform in the bulk; 

h(-,s) = /(s) in O. (1.34) 

As a consequence, the concentration held satishes the following differential inclusion: 

/(s) G /So sign {dsu) + ku, (1.35) 

which is known to exhibit hysteresis (actually, it reproduces the well-known stop operator, 
see e.g. pi25]). 

The initial boundary value problem. We consider the evolution in a smooth domain 
in the expanse of time (0,T). We suppose that at time t = 0 the concentration held 
be given by a prescribed function Uq{x), x E Q. We also suppose that a time-dependent 

chemical potential p\,{x,t) be prescribed on for all x G F at all times t G [0,T]. 

At each particular time, we harmonically extend p\) to the interior of (still denoting 
by p\, the harmonic extension) and we introduce the characteristic time and lengthscale 
To = q; and Lq = hf0 as well as the new variables and functions: 

t = x = ^, p = p-p^, ^(r) =/3(r/To). (1.36) 

Lq -^0 

^We assume that all energies densities per unit volume are renormalized to a reference value, and so 
are dimensionless. 













Regularization of a forward-backward parabolic equation 


We express the system in terms of these new variables and we henceforth drop tildas, so 
as to obtain the following problem: 


dtu = A/i, I 

p = dtu + (u) > inf2x(0,T) (1.37a) 

^ e (dldtu) ) 

with the boundary condition 

p = 0 on r X (0,T), (1.37b) 

and the initial condition 

m(-,0) = 'Uo in fl. (1.37c) 

Before proceeding, let us set a precise outline of the paper. In the next section, we 
introduce notation, assumptions on the data of the problem and state the main existence 
and uniqueness result, which is complemented by the continuous dependence of the solu¬ 
tion with respect to the data of the problem. In Section [3] we proceed by exploiting the 
a priori estimates on the solutions of the system we use to prove existence and regularity. 
In Section m we show the continuous dependence estimates on solutions. Finally, in Sec¬ 
tional we provide a detailed proof of the existence and the uniqueness of the solution. To 
this aim, we hrst establish our result for a “regularized” version of the system, obtained 
by a suitable truncation of the coarse-grain free energy mapping; then, owing to the esti¬ 
mates carried out in Section [3l and by making use of a maximum-principle argument, we 
establish our result for the original problem. 


2 Notation, assumptions, and main results 

Before stating the problem we are dealing with and the main existence result, let us 
make precise our assumptious ou the data aud the notatiou we use. In the sequel, 17 is a 
bounded smooth domain in with smooth boundary F. We introduce the spaces 

H := L^(17), V := i7o(17), W := n Ro(17). 

We endow H, V, and W with their usual scalar products aud norms, and use a self- 
explanatory notation, like || • ||y. For the sake of simplicity, the same symbol will be used 
both for a space and for any power of it. We note that the norms ||u||y and ||Vu||ir are 
equivalent for v G V, thanks to the Poincare inequality. In addition, let us point out 
that, after identifying H with its dual, the triplet (V, H, V) is a Hilbert triplet (where 
V coincides with the Sobolev space Hence, we use the notation (•, •) for the 

duality pairing between V and V. Given a hnal time T > 0, we set 

Q:=Qx(0,T), S:=Fx(0,T). 

As far as the data of the problem are concerned, we assume that the domain of ^|J is an 
open interval (a, b) C M, where a and b could be taken equal to —oo and -|-cxo, respectively. 
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and require that 


e C^(a, b), 

> 0 for all r G (a, 6), 
lim 'ijj'{r) = — oo, lim 'ijj'{r) = +oo, 

r—)-a+ r^b~ 

%lj"{r) > —Ki for all r G (a, 6), 


for some positive constant Ki, note that in fl 2 . 1 cp has to become —oo if a = - 
b~ reduces to +oo if & = +oo.For the initial concentration uq we suppose that 

Uq E H, 3ao > a, bo < b such that Oq < Uo{x) < bo for a.a. x E Q, 

whence both uo and 'ip'{uo) lie in L°°{Q). Concerning the known datum pt,, we 
that 

Finally, as to the (possibly) multi-valued mapping /3, we let (see [3]) 

(3 = dC, with ^ : M —)■ [0, -|-oo] convex and 
lower-semicontinuous, such that C( 0 ) = 0 . 


(2.1a) 

( 2 . 1 b) 

( 2 . 1 c) 

(2.1d) 

-OO and 

(2.1e) 

assume 

(2. If) 


( 2 . 1 g) 


Let us specify a weak formulation of the problem in the set of the Hilbert triplet 
(C, H, V). Let us dehne the operator A -.V ^ V, corresponding to the “weak realization” 
of the Laplace operator —A (combined with homogeneous Dirichlet boundary condition) 
in the duality between V and V, by letting 

{Avi,V2) ■■= / Vvi-Vv2dx, Vi,V2EV, ( 2 . 2 ) 

Jn 

and specify its inverse A~^ : V' ^ V, such that for w,z E V' and v E V 

{Av,A~^w) = {w,v), {w,A~^z) = {z,A~^w) = f V(H“^t(;) ■V{A~^z). (2.3) 

Jn 

We introduce a norm in C', denoted by || • ||*, which is equivalent to the usual one, 

\\w\\l = {w, A~^w), weV. (2.4) 

Definition of solntions. We say that a triplet (m, p) is a weak solution to the problem 
flL^ if 

u E C'([0,T];iL), e G C°([0, T]; if), /i G C°([0,T];C), (2.5a) 

ij\u) EH\d,T-H) (2.5b) 

and the following equations are satished: 

dtu{t) + Afv(t) = 0 in W, for all t E [0,T], (2.6a) 

fi{t) = {dtu + ^ + fi\, +'ij^'{u)){t) a.e. in H, for alH G [0,T], (2.6b) 

^{t) E (3{dtu{t)) a.e. in H, for all t E [0,T], (2.6c) 

m(0 ) = Uo a.e. in H. ( 2 . 6 d) 
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Regularization of a forward-backward parabolic equation 


Theorem 1 (Existence and uniqueness). Under the assumptions fl2.ip . there exists a 
unique solution to the problem fll.371) . in the sense we have speeified above. In addition, 
it results that 

u e L^{Q) (2.7) 

and the solution is strong, that is, ja G C°([0, T]; W) and equation fl2.6al) ean be replaced by 

dtu{t) — Afi{t) = 0 a.e. in Q, for all t & [0,T]. (2.8) 

Moreover, a continuous dependence on the data holds: namely, if p.i), {u 2 ,^ 2 ,h^ 2 ) 

are two solution triplets corresponding to the initial data uqi, Uq 2 and bulk data fi\,i, ia\, 2 , 
respectively, then their difference satisfies 

11^1 — it2|lci([o,T];R) + ll'^i “ 'C2||co([o,t];R) + WUi ~ U2\\cO([0,T]-,W) 

< R (ll/^bl ~ U\)2\\cO{[0,T]-,H) + 11^01 — W02||r') (2.9) 

for some constant R depending only on the structural assumptions stated in (Bill. 


3 Basic estimates 


In this section, for reader’s convenience, before proving Theorem (H we recover the a 
priori estimates we can derive on the solutions of fl2.6p . Note that we are using the same 
notation C for possible different positive constants depending only on the data of the 
problem. 


Energy estimate. We hrst show that system 02.61) admits in a natural way a so-called 
“energy estimate”. Indeed, once u and /i are solution components for the problem 01.37p . 
we are allowed to test 02.6al) by p and 02.6bD by dtu. Then, we combine the resulting 
equations and integrate by parts in time over (0,t); by exploiting Holder’s and Young’s 
inequalities and using the chain rule and the smoothness of we hnd that 


/ il){u{t))dix+ / f dtU + IVdxds 

Jn JoJn 

< / fi{uo)dx+ / / \fi\,dtu\dxds, 

Jn JoJn 

< f fi{uo)dx + ^ f f l/Xbpdxds-F^ f f \dtu\‘^dxds. 

Jn Jo Jn 2 Jq Jq 


(3,1) 


Owing to the monotonicity of jS and the fact that 0 G /3(0), we have that 


f dtU dxds > 0. 


0 Jn 


(3,2) 


As 02.5ap holds, we also point out that u{t) = uq + dtu{s)ds, whence 


\nml < nuofn + 2T 


||atu(s)||^ds. 
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thanks to the Holder inequality. Thus, by virtue of the Poincare inequality and the 
nonnegativity of ^|J as well, fl3.ip and fl2.ip yield the estimate 

||'0('*^)IU°°(O,T;Ll(O)) + \\u\\m(o^T-,H) + ||hlU2(0,T;y) < C. (3.3) 

Note that fl3.3p entails in particular that u lies between a and b almost everywhere in 
Q.Hence, recalling that and P are smooth enough and that p has null trace on P, by 
a comparison in fl 2 . 6 al) and standard elliptic regularity estimates we obtain, in addition, 
02.81) and 

11^11^2(0,T;W) < C. (3.4) 


estimate for chemical potential. Let us hrst introduce po G H as the unique 
solution of the nonlinear elliptic problem 

+ (-^ +/^) ^(ho ~ hb(0) —'0 ^(mo)) = 0. (3.5) 

Indeed, note that /it>(0) + 'ip'iuo) E H hy 02.1fp . 02 . 1 ep and 02 . 1 al) : then, concerning the 
sum of the two maximal monotone operators A (when restricted to W with values in 
H) and 

n i-A (/ + - /ib(0) - 'i/j'iuo)) 

we can apply [3l Cor. 1.3, p. 48], which ensures that the sum is maximal monotone and sur¬ 
jective thanks to the Lipschitz continuity of the second operator and the coerciveness of A. 
Moreover, the uniqueness of po G IT solving 03.51) follows from the strong monotonicity 
of A. Clearly, from 02.5p and 02. 6 p we have that /i(0) = po and 

dtu{0) = Apo =: Uq E H and ^(0) = po - Mq - hb(0) - i>'{uo) =: G H. (3.6) 


We emphasize that satishes G (3{ uq ) almost everywhere in Next, in order to 

show the L°° estimate for /i, let us £x t G (0,T] and, for n G N, set = t/n, = ir„, 
i = 0,1,..., n. For typographical convenience, we henceforth omit the dependence of t" 
on n. In view of 02.5p . equations 02.6ap and 02.6bp hold at the times ty. 


dtuiti) H/i(C) = 0 in W, 

/i(C) = dtu{ti) + ^{ti) + (/ib + 'ip'{u)){ti) in H, 


(3.7) 

(3.8) 


for i = 1,..., n. Test 03.7p by p(ti) — p(fi_i) and the difference of equalities 03.8p at the 
steps i and i — 1 by dtu{ti). Then, by combining the results is not difficult to check that 

V K/ K/ 

+ \\\SMU)\\1 - + \\\dMti) - dMU-M 

+ (^(C) - ^(C-i), dtu{ti)) + + 'ijj'{u)){ti) - (pb + dtu{ti)) = 0, (3.9) 

where (•, •) denotes the scalar product in H. Now, by the properties of the sub differential 
(3 = d( (which is a maximal monotone graph) it turns out that the inclusion (see 02 .6cl) 
and ( |2.1g|) ) 


^(C) G dC{dtu{ti)) 
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can be rewritten as 

dtuiti) e dC{i{ti)) (3.10) 

almost everywhere in where •= sup^gR(nt<; — ({v)), tc G M, is the Legendre- 

Fenchel transform of (, and its sub differential d(* coincides with the inverse graph 
of (3. Then, using the dehnition of subdifferential, it is straightforward to infer that 




Hence, on performing summation in fId.Op for i = 1,..., n, we plainly deduce that 
^ Kt)\^dx + ^\\dtu{t)\\jj + j C*(e(t))dx 

k) k) 

dtu{ti)) 


2 = 1 

1 


< » / |Vw|Mi + -||ui||i+ / C({o)da:. 


(3.11) 


Owing to 112.1 til and Il 2 .. 3 hll . it is a standard matter to infer that 

^ + nn)m - (it, + it'(a))(t.-.) ^^ j-r + 

Ti Jo Jn 

as n —)■ cxo, and assumptions (12.1 dp and fl2.1fp easily yield 

^t(hb +'Ip'{u))dtU dxds > —\\dtlJ.\,\\L 2 { 0 ^T-,H)\\dtU\\L^{ 0 ,T-,H) — -^l||^t'w|lL 2 ( 0 _T;n) — 



0 JQ 


the last inequality being due to the previous estimate fl3.3p . As .^o ^ dCiu'o) and conse¬ 
quently C*('Co) + C('*^o) = 'Co'i^o almost everywhere in H, from p.lgp it follows that 


C*(^o)dx < ll^ollHllMollfi - / C(wo)dx<a 


(3.12) 


On the other end, it is easy to check that (* is non-negative, whence 


/ r(e(t))da; > 0 . 

Jn 

Then, passing to the limit as n —)■ oo in fl3.1ip and exploiting the previous remarks, we 
hnd out that 

\\Kt)\\l + \\dMt)\\l<C for all t G [0,T]. (3.13) 

At this point, by comparison in fl2.6ap and thanks to well-known elliptic regularity results 
combined with the Sobolev embedding W C L°°{Q), we recover 

||hllL°°(0,T;W) + l|hlU°°(Q) < C. (3.14) 

In particular, because of the boundedness of /it, postulated in Assumption fl2.1fp . we 
infer that there exists a constant M such that 


111*-< M. 


(3.15) 
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L°° estimate for concentration, 
problem in H: 


We combine fl2.6bj) - (l2.6dj) 


to obtain the Cauchy 


dtu{t) = {I + P) - ixPt) - ijj'iuit))), te[0,T], 

m(0) = Uq. 


(3.16a) 

(3.16b) 


Now, by assumptions (12.1 cp and fl2.1ep there exists two constants k^,k* e (a, b) such that 


ip'{r) > M for all r > k*, (3.17a) 

'p'{r) < —M for all r < fc*, (3.17b) 

k^ < ttQ < Uq{x) <bo<k* for a.a. a: G 11. (3.17c) 

We test fl3.16ap by {u — k*)^ and —{u — k^)~, then integrate over (0,t). We note that 
(/ + P)~^{r) has the same sign of r and the right hand side of fl3.16ap is nonpositive if 
u > k* and nonnegative if m < A;*, thanks to fl3.17p and fl3.15p . Then, after integration by 
parts in time, it is a standard matter to infer that 

k^ < u < k* a.e. in Q, (3.18) 


which entails (1^ . 


4 Continuous dependence on the data 

Consider a pair of data {woi, /ibi}; * = 1)2, fulfilling (I2.1ep . (I2.1fp and let 


i = 1,2, be the corresponding solutions. Then, the triplet with u := Ui — U 2 , 

^ ■= — ^ 2 , ■= k'l — /^ 2 ) satisfies (cf. (12.6p and (I2.8P ') 

dtu{t) — Afl{t) = 0 a.e. in 11, for all t G [0,T], (4.1a) 

fi{t) = {dfU + ^ + fib +— 'ijj'{u2)){t) a.e. in H, for all! G [0, T], (4.1b) 

^i{t) G P{dtUi{t)) a.e. in H, for all t G [0,T], i = 1,2, (4.1c) 

11(0) = uq a.e. in H, (4-ld) 


where uq = uqi - mo2 and /ib = /^bi - hb2- 

In view of the regularities in (I2.5ap . we can test (14.lap by (I4.1bp by dtu(t) and 
add the resulting equations. In particular, by virtue of (I4.1cp and by the monotonicity of 
P, we have that f^^(t)dtu(t) > 0 and consequently 

[ |V/i(l)pda;+ [ dx < f {\fiPt)\ + - 'ilj'{u2{t))\)\dtu{t)\dx. (4.2) 

Jq Jq Jq 

Now, since ^|J is twice continuously differentiable, its derivative is locally Lipschitz- 
continuous. Moreover, by the estimate (I3.18p . it turns out that both solutions Ui stay in 
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Regularization of a forward-backward parabolic equation 


a bounded interval J. Consequently, we have that — '0X^2 )! < ||'0^1Il°°(j) 1^1 and, 

by Young’s and Poincare’s inequalities, we infer 

mmi + \\ 9 Mmi < c(iirt(()iii,+ii«wiii,). (4.3) 

Given that ||h(t)|||^ < C ^||ho|||^ + /q ||9ih(s)||f^ dsj, from fl4.3p it follows that 

mmi + mmi 

< c (^\\Mt)\\H + \\Mt)\\H + ||aih(s)||^ds^ for alH e [0,T]. 

Thus, an application of the Gronwall-Bellmann inequality yields 

|lh|lcO([0,T];U) + \\dtU\\cO{[0,T]-,H) < C (||/lb||c'0([0,T];_ff) + llholln) • (4.4) 

Then, the analogous estimates for Ajl (and consequently for jl in G°([0,T]; W)) and ^ in 
C^{[0,T]]H) follow from fl4.4p by a comparison in fl4.1ap and M.lbli . which helps us to 
conclude the proof of fl2.9p . Of course, 02.91) implies in particular the uniqueness of the 
solution to the problem 02 .6p . 


5 Existence of solutions 

In this section, we give details on the proof of the existence of the solution to our problem. 
We use a contracting argument. First let us make a truncation of the function -0, which 
allows us to exploit the above a priori bounds on the solutions to fl2.6p . Let K and k* 
be two constants fulhlling 03.171) . It is not hard to check that, thanks to the assumption 
02.1cl) . there exist constants and K* such that 





(5.1) 

and 

<(i^*) > 0, > 0. 

We introduce the following truncation of ip: 


(5.2) 

ip*{r) = 1 

{'tppr) 

pjiK*) 'ip'{K*){r - K*) l'ip'\K*){r - K*f 

yiK^) + 'iP'{K^){r - IQ + iniQir - KQ 

if K^<r< K\ 
if r > K*, 
if r < JF*, 

(6.3) 

and we denote by 

L ;= max '0*(r) < -|-cxo 

rGK 


(5.4) 


the Lipschitz constant of its derivative -01. We note on passing that the truncated function 
-0* satishes the assumptions 02.lap - 02.Idp with (a, 6) = (—cxo,-|-cxo). In particular, the 
bound from below02 .1dp holds for '0" with the same constant —Ki as for ip". 


Next, we consider the set 

S-- {K6C"(|0,T|;i/): t>(0) = «„} 


(5.5) 
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and we introduce the map : S ^ S which to every v E S associates u = J^{y) dehned 
by 

u{t) = J^{v){t) := Uo+ [ (/ + /3)"^(/i(s)-/ib(s)-'0*(w(s)))ds, te[0,T], (5.6) 

Jo 

where fi(t) denotes the unique element of V that solves the nonlinear elliptic equation 

ylp(t) + (J + /3)"^(/i(t)-/ib(t)= 0, te[0,T]. (5.7) 

Before proceeding, let us comment on the existence of a unique p G C°([0,T];1/) 
satisfying 05.711 for some v hxed in S. First, we recall 02.1fjl i and observe that is 
Lipschitz continuous, so that the function t i-A- is continuous from [0,T] 

to H. Then, for all t G [0,T] there exists a unique p(t) fulhlling 05.7p : this can be shown 
arguing as for 03.511 and using [31 Cor. 1.3, p. 48]. Moreover, as (/ + /5)“^ is monotone and 
Lipschitz continuous, it is not difficult to check that fi G C°([0, Tj; C). Once jj, is found, 
the function u = JF{v) E S is completely determined from 05.6p . 

Eventually, our aim consists in applying a hxed point argument: indeed, we will see 
that any hxed point for the operator JF turns out to be a solution to the problem made 
precise by 02.5^ - 02.611 . To this aim, we are going to show that some power JF'j (j G N) is 
a contraction mapping in S and, as a consequence, it admits a unique hxed point, which 
results at the end to provide the unique solution to our system. 

To this aim, we pick a pair {ni}i=i ,2 C S and set 

Ui := := Pi - pb - - dtUi, 

where p* is the solution to 05.7p corresponding to Uj, f = 1,2. Then, it is easy to verify 
that 

dtUi{t) + Afii{t) = 0 in V', for all t E [0,T], (5.8a) 

^i{t) = (dtUi + + ^jJ'{vi)){t) a.e. in for alH G [0, T],, (5.8b) 

^i{t) E (3{dtUi{t)) a.e. in for all t E [0,T], (5.8c) 

Mi(0) = uo a.e. in (5.8d) 

for i = 1,2. Now, we use the notation u for the diherence of ui — M 2 , and the same 
notation for jl and v. We take the diherence of 05.8all for i = 1, 2, test it by A~^{dtu(t)) 
and, at the same time, we test the diherence of 05.8bp by dtu{t). Then we combine the 
obtained equalities and use the properties of A~^ stated in 02.311 and 02.411 . Hence, we 
have that 

\\dtu{t)\\l + WdMtWn + [ ^{t)dtu{t)< [ -'i/ji{v2{t))\\dtu{t)\ 

Jn Jn 

< hlMml + \l W.Mt)) - ’P'.i'Mi))?- (5.9) 

k) 

Due to 05.4p . we handle the right hand side observing that 

Jn Jn 


(5.10) 
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In addition, by the monotonicity of (3 and fl5.8cp we dednce that 

[ ^{t)dtu{t) > 0 . 

Jn 

Thus, we easily obtain 

\\dtu{t)\\H < L\\v{t)\\H (5.11) 

and, as both ui and U 2 satisfy the same initial condition fl5.8dD . we can easily infer that 

\\J^{vi{t)) - T{v2{t))\\H < L [ ||ni(s) - n 2 (s)||Hds for all f e [0,T], (5.12) 

Jo 

This inequality leads to 

||J'(ni) --T(n 2 )||c 0 ([ 0 ,t];R) < - ^' 2 ||co([o,t];R) for all t e [0,T]. (5.13) 

An iteraction of the argument, due to 05.121) and 05.13p . leads to 


\T\v,)-J^‘^iv2)\\coim,H) <L ||^(ni)(s) -^(n2)(s)||/^d5 

. LtT 


s\\vi - T2||c'0([0,5];R)ds < “ ^2 ||cO([ 0 ,t];R) • 


'0 


By iterating j times, we find || J'^(ni) - {v 2 )\\co{[o,t]-,H) < ^\\vi - V 2 \\co{[o,t]-,H) for all 

t G [0,T], whence, in particular. 


- J^\v2)\\cO(10,T]-,H) < 


(LTy 

j! 


11^1 — U2||cO([0,T];n)- 


(5.14) 


Thus, for j large enough turns out to be a contraction mapping from S into itself, as 
announced; as a consequence, has a unique fixed point m*, which is also the unique 
hxed point for T. In view of 05.81) . this hxed point yields the triplet {u*, ^*, ^*) that solves 
the problem 02.61) in which ijj is substituted by '0*. Of course, for {u*, fi*) we can repeat 
the estimates carried out in Section [31 In particular — and this is the crucial point — we 
can derive for fi* the same estimate as 03.15p . namely. 


ll/i* - /i[,|| 2 ,oc(Q) < M, 


(5.15) 


with the same value of the constant M. In fact, if one checks carefully the estimates, one 
can see that appears in 03.11) with the integral of 'ip*{uo) = ^Jiuo) and in 03.6p with 
'^((mo) = fj'iuo) (cf. 02. lei) and 03.17cp h and especially with the constant Ki in 02.1dp 
which, as observed at the beginning of this section, can be the same for and V’*- In 
addition, by its very dehnition the derivative satisfies 


-^((r) > M for all r > k*, 
V'Ul) < —M for all r < k^, 


(5.16a) 

(5.16b) 


as xf' does in 03.17ap - 03.17bD . Thus, a repetition of the argument leading to 03.181) yields 

k^ <u* < k* a.e. in Q (5.17) 
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and, since -0 = V’* k*]) il5.17p entails 


a.e. in Q. 


(5.18) 


In other words, is actually a solution to the original problem fl2.6p and ful- 

hlls fl2.7p . Moreover, is the unique solution of fl2.6p . owing to the continuous 

dependence property fl2.9p proved in Section 01 Finally, recalling the smoothness of 12 and 
r and the homogeneous boundary condition on F, by fl2.5p . fl2.6bp and standard elliptic 
regularity estimates we obtain fl2.8l) and the regularity /i G C°([0,T];hF) for 
Therefore, Theorem [T] is completely proved. 

Remark 1. Note that if we assume Uq G H^{Q) and /i), G L^(0, T; 77^(12)) besides fl2.1el) 
and fl2.1fp . then we can recover the additional regularity u G 77^(0, T; 77^(12)) for the 
solution component u. Indeed, it suffices to take formally the gradient of fl3.16ap and test 
it by V{dtu). Thanks to the Lipschitz continuity (with Lipschitz constant 1) of (7 + 
and of ipl (which can replace as we have seen) with constant L (cf. fl5.4p h by the Young 
inequality we easily infer that 


1||v(a,«)(«)f„ < ||v(fi- p.Mt)\\l + 


Now, pointing out that E L^(0, T; 77^(12)) (cf., e.g., fl2.6ap i. as 



we can easily apply the Gronwall lemma and hnd out that dtu G L^(0, T; 77^(12)). 
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